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["T I ■ distance at galaxies is then evaluated by using the SBF method. For this pur- 

00 ■ pose, theoretical SBF magnitudes for the ACS@HST filters are computed for 
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us obtain i/o = 76 ± 6 (statistical) ±5 (systematic) km/s/Mpc. 



distance scale 



1. Introduction 

The Surface Brightness Fluctuation (SBF) method is a powerful technique to derive dis- 
tance to galaxies when individual stars cannot be resolved. The method was first introduced 
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by iTonry fc Schneiderl (119881 ) and is based on a simple concept: the poissonian distribution 
of unresolved stars in a galaxy produces brightness fluctuations between pixels of the galaxy 
image. The resulting pixel-to-pixel variance of the fluctuation is inversely proportional to the 
square of the distance, and therefore can be used as distance indicator. Although limited by 
pixel resolution and high S/N ratio, the SBF method showed the capabil ity of spanning an 



interval of distances extremely wide, ran ging from Local Group obje cts (lAihar et al 



Raimondo et al.ll2005l : iRekola et al.ll2005l ) up to galaxies at 100 Mpc ( jThomsen et al. 
Jensen et al.l 2001 ). 



1994 



1997 



Ideally, SBF measurements can be derived for almost any morphological type of galaxy, 
provided that the region analyzed has a nearly regular /smooth luminosity proflle. Indeed, 
the presence of morphological irregularities represents a constrain to reliable measurements 
of SBF, for example, if dusty regions are present, the application of SBF technique requires 
an accurate masking of dust. 



As recognized by iTonry et al.l (Il990l ). the properties of the stellar populations in the 
galaxy cannot be neglected in calibrating the absolute SBF magnitude. This is usually 
overcome by providing a relation between the SBF magnitude and a galaxy color. Once such 



lenticular and sp iral galaxies with prominent bulge (e.g. 



Tonrv et al. 


2001; 


Mieske et al.l 


2006 



Mei et al.l 120071 ). Moreover, the SBF technique has b ee n successfully 



brightness dwarf ellipticals (iJerien et al. 
Mage llanic Clouds stellar clusters (lAjhar et al. 
2005h . 



1998, 2000a 



1994 



b; 



Rekola et al 



Gonzalez et al 



appl ied to low surface 



2005), and Galactic or 



2004; Raimondo et al. 



An interesting case is represent ed by shell el l iptica l galaxies, a small number of which is 
included in the sample observed by iTonry et al.l (120011 ). Shell structures are generally faint 
and sharp stellar f eatures, and they are considered a robust in dicator of past merging or 
interaction events (IMalin fc Carterlll983l ; IWilkinson et al.ll2000l ). The stellar population in 
the shell depends on the galaxy with which the rnerger has taken place a nd on the time spent 



since the shell structure has formed (lPencelll986l ; IWilkinson et al.ll2000 ) . Many authors have 



found shells redder or blu er than the underlying galaxy (IForbes et al 
I999I : ISikkema et al.l[2007h . 



1995; TurnbuU et al. 



It is reasonable to expect that the presence of shells might influence the SBF signal of 
the galaxy. In this respect, the high quality of Advanced Camera for Surveys (ACS) images 
is crucial to remove the shell structure, allowing the measurement of SBF of the galaxy, even 
at high distances. Thus, the SBF measurement remains a relevant tool to investigate the 
distance of shell galaxies. 
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In the present work, we derive SBF magnitudes of a sample of distant shell-ellipticals 
taking advantage of the high resolution power of the ACS on board of the Hubble Space 
Telescope (HST). We select deep ACS images and measure F814W-band (~ I-band) SBF 
of the four galaxies: PGC6510, PGC 10922, PGC 42871 and PGC6240, with recession ve- 
locities reaching 8000 km/s. The main goal of our work is to measure their distance by 
using the SBF technique. No previous distance determination is available for any of these 
galaxies, except for the kinematic distance modulus. The farthest object of the sample, at 
a distance of ~ 100 Mpc, lies near the upper limit o f distances obtained using optical SBF 
measurements (IThomsen et al.lll997l : iMei et al.ll2003l ). It is worth noticing that none of these 
objects has been previo usly analyzed in detail, with the only exception of the recent work of 
Maybhate et al.l (120071 ) on the globular cluster system of PGC 6240. 



In the next section we present the observational data and the general properties of the 
selected galaxies. The procedure to derive the SBF measurements and the data analysis 
are described in Sect. [31 The results are presented and discussed in Sect. HI where we 
provide new SBF predictions together with new calibrations of absolute SBF magnitudes for 
the ACS@HST filters based on single-burst stellar population (SSP) models. In the same 
section, an evaluation of Hq is also obtained. A brief summary is given in Sect. O 



2. The observational data 

The HST images used in the present work were obtained with ACS in its Wide Field 
Channel mode. The large field of view, high resolution and sharp point-spread function (PSF) 
characterizing ACS images are critical requirements to attempt the SBF measurements of 
the distant galaxies we have selected. 

The images were retrieved from HST archive. A requirement in selecting the data is 
that the SBF S/N ratio in the F814W band is ^ 5, as suggested bv lBlakeslee et all Jl999[ ). 



This condition is verified for the three ellipticals PGC 6510, PGC 10922, PGC 42871, while 
PGC 6240 images have slightly lower S/N (~4). For these galaxies F475W images were 
also available, though the exposure times prevented their use to measure SBF (S/N << 
5). However, they were retrieved and analyzed to obtain integrated magnitudes and colors. 
The images are associated with proposal GO 10227 (PI: P. Goudfrooij) designed to study 
the globular cluster system of the four giants, post- star burst shell ellipticals. The main 
properties of the selected galaxies are reported in Table [H 

We downloaded the ACS images processed with the standard calibration pipeline (CALACS), 
that includes bias, dark and fiat-fielding corrections. The images from the archive still re- 
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quire the final image combination and the correction to remove the geonietric distortions. 
To this purpose we used the PYRAF task multidrizzle flKoekemoer et aLll2002l ). that also 
provides the bad pixel identification and the cosmic-ray rejection. No sky subtraction was 
performed at this stage. To reduce the effect of noise correlation introduced by the drizzling 
procedure (IFruchter &: Hookl 1200211. we adopted the L ANCZ0S3 kernel which, as demon- 
strated by Mei et al. ( 2005 ) and ICantiello et aL ( 2005 ). is adequate for the purpose of SBF 
estimation. 



3. Data Analysis 



T o deriv e the SBF measurements we follow the same basic procedure adopted in lCantiello et al 



(120051 . 120071 ). In this section, we summarize the relevant steps to measure SBF magnitudes 
and we enlighten the differences from the quoted works. 

First, a provisional sky value is derived from the corner with lowest sky counts, and 
a mask of the bright sources (saturated stars, extended galaxies, etc.) is obtained. Then, 
the brightness profi le of the galaxy is modeled using the IRAFO/STSDAS task ELLIPSE 
(jjedrzejewskil 119871 ). After the galaxy model is subtracted, we derive a new mask of the 
faint sources which clearly appears after the galaxy light is removed. This new mask is fed 
to ELLIPSE to obtain a new galaxy model. 

We use the isophotal geometry, as derived with the last galaxy model, to get the surface 
brightness and color profiles of the galaxy. The galaxy profile is fit using a de Vaucouleurs 
law to find the sky as the zeropoint constant in the fit. For all galaxies, we find that a de 
Vaucouleurs law is well suited in the 1-5" region of the galaxy. This result confirms the 
known eviden ce that the photomet ry of shell galaxies is the same expected for a normal 



elliptical (e.g. IWilkinson et al.ll2000l ). 



The new sky value is adopted and all the above steps are iterated, till convergence. At 
the end of these iterations we have the final sky value, the mask of external sources, and the 
best galaxy model. The final surface brightness and color profiles for the four galaxies are 
shown in Fig. [H It is worth noting that the average color of the four galaxies is bluer than 
in normal ellipticals, and the color gradient is po sitive or nearly flat. Both character istics 
can be considered normal for this class of galaxies (ITamura et al.ll2000l : iLee et al.ll2006l ). For 



"^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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all da ta we apply a K-correction Kpsuw ^ Kj ^ 0.5 ■ z and -ft'F475iy ^ Kg 2 ■ z (iPoggianti 
19971 ) for integrated co l or. T he extinction correction are evaluated using the prescriptions 
given by lSirianni et al.l (120051 . their Table 14). 



We subtract the sky value and the galaxy model from the original image. This operation 
leaves a large scale residual background, due to mismatch of the real galaxy with the model. 
The large scale residuals a re removed using the ba ckground map derived with the photom- 
etry package SExtractor (IBertin &: ArnoutsI 119961 ). We have carried out some numerical 
experiments to determine the best background map parameters able to provide both a good 
subtraction of the large scale residuals, and the best possible removal of the shell features. 
After several experiments we have chosen a mesh size of 15x15 pixels^ (BACK_SIZE=15), 
with 3 background-filtering meshes (BACK_FILTERSIZE=3). 

Up to this point, the procedure is applied to both filters. The sky, galaxy model and large 
scale residuals are subtracted to the original frame to derive the residual frame. Note that 
possible dust regions, as found in PGC 42871, are masked out since they could compromise 
the SEP measurement. Such regions are better recognized from the P475W image, the same 
dust mask is used for both filters. 

We run SExtractor on the residual frame to obtain the photometric catalog of external 
sources (foreground stars, globular clusters and background galaxies). This photometric 
catalog is used to construct the luminosity functions (LPs), which will be used to estimate 
the contribution to the fluctuations coming from faint unmasked external sources (P^es)- 



The procedures used to fit the LPs are the same described in ICantiello et al.l (120051 ). 
In particula r, we adopte d M/ = —8.5 mag as the absolute Turn Over Magnitude (TOM) of 
the GCLP (lHarrisll2001l). while th e exponent for the power law LP of background galaxies 
is 7 = 0.34 (jBernstein et al. 2002). These parameters are used to start an iterative fitting 
process where the surface density of galaxies and globular clusters, and the galaxy distance 
are allowed to vary until their best values are found via a maximum likelihood method. Pig. 
|2] shows the best fit of the observed LPs obtained for the four galaxies, used to derive Pres ■ 

Before computing the power spectra of the residual image and of the PSP, and evaluate 
their azimuthal averages, we mask regions with residual contamination from shells, and 
divide the residual frame by the square root of the galaxy model. It is worth noting here 
that only the shells not completely removed by the background map subtraction, i.e. the 
most prominent shells, required further masking. 

The frames obtained with the analysis described up to this point are shown in Pig. [31 
We point out that in general the residuals appear similar to typical residual frames of SBP 
measurements of normal ellipticals. Nevertheless, in Pig. [3] we adopted a scale to emphasize 
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the regions where the procedure described failed in the complete subtraction of shells (in 
particular, PGC42871 and PGC6240). To quantify the additional systematic uncertainty 
due to unsubtracted shells, we performed a specific test described at the end of this section. 

Next, the total fluctuation amplitude is determined as the constant factor Pq multiplied 
by the PSF power spectrum, to match the power spectrum of the residual frame P(k): 



Pik) = Po-E{k) + Pi, 



(1) 



where Pi is the white-noise component, and E{k) is the convolution between the PSF power 
spectrum and the power spectrum of the mask function. The latter mask also takes into 
account the shape of the galaxy annulus adopted for the SBF measurement. The minimum 
annular radius has been fixed to be the minimum radius without dust contamination, while 
the maximum radiu. i= fixed to the region where the galaxy to sky counts ratio is > ig 

A robust linear least-squares method JPress et al.l Il992h is used to fit Eq. ([T]). The 
lowest k-numbers j k < 250), that have b een corrupted by the subtraction of the smooth 
background profile (IBlakeslee et al.l 119991 ). and the high k-numbers {k > 600), that have 



been corrup ted by the drizzling procedure, are excluded from the fitting (see ICantiello et al. 
2OO5I . I2OO7I . for more details). 



Fig. m exhibits the azimuthal average of the power spectrum for each galaxy, the best 
fit lines are also shown. 



Finally, the SBF magnitude is evaluated as: 

mpsuw = -2.5 log (Po - Pres) + m* - Apsuw - K 



F814VK 



(2) 



where Pres is the extra contribution of unmasked external sources, evaluated from the fitted 
LF as described in the quoted papers. The amplitude of Pres is small for all galaxies of the 
sample, being on average Pres/Po ^ 0.08 (Table [2D- m* is the zero-point ACS magnitude in 
the VEGAMAG system reported by ISirianni et al.l (l2005l . rrfp^^^^-yY = 25.501 mag), Apsuw 
is the extinction correction in the F814W passband, and Kp^^^^r is the i^-correction term. 
We apply Kp^^^y^ ^ 7 ■ z after iThomsen et al.l (119971 ) . 



^ As ICantiello et al.l (|2007l ). we consider one single annulus per galaxy, because of the lower S/N ratio of 
the images, and because of the, on average, small available area. It must be noted that in our procedure the 
contribution of the external sources is evaluated taking into account their radial position, so that using one 
single annulus does not introduce any systematics on Pres- Moreover, we have carried out a three-annuli 
measurement on the galaxy with the largest spatial extension, PGC6510, as a result we found that the 
averaged value agrees within uncertainties with the measurement reported in Table [H 
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We used ISirianni et al.l (120051 ) equations to transform the F814W and F475W into the 
standard B and I magnitudes. However, in the forthcoming section we will also take into 
account the magnitudes in the ACS photometric system. 



T he ma in differences between the present data and the images used by ICantiello et al. 



(I2OO5I . I2OO7I ) are i) the much greater distance of the objects, which affects the LF fitting, 
and a) the rather complex shell structure. For this reason we have performed two additional 
tests with respect to the ones described in the quoted papers. 

The first one is related to the Globular Cluster LF (GCLF). The data quality of the 
images in some cases does not allow to reliably sample the T OM of the GCLF. This is 
confirmed by the recent work published by lMaybhate et al. J2007h on the GCLF of PGC 6240, 
based on the same data used here. These authors do not provide an estimation of the TOM, 
however, from a visual inspection of their Fig. 11, the extrapolated TOM appears ~1 
magnitudes fainter than the one adopted here. Since the large uncer tainty, we perfor m a 
specific test to evaluate its effect on the SBF. Adopting the GCLF of iMaybhate et al.l . the 
SBF magnitude changes a few hundredths of mag ( ^ 0.05 mag), that is a ~ 3% uncertainty 
on the distance. 

The second test concerns how the presence of shells affects SBF. To evaluate this effect 
we selected the farthest galaxy in the sample (PGC6240), then we run the whole SBF 
measurement procedure on the image, with and without masking the most prominent shell. 
The result of this test shows that the final SBF magnitude changes by ^0.25 mag. Even 
though this is not an exhaustive evaluation of the uncertainty caused by shell features, it is 
a robust suggestion that the adopted procedure should keep this source of uncertainty fairly 
below 0.25 mag. Thus, we consider this value as the upper limit of systematic uncertainties 
due to possible unmasked residual shells. 

In conclusion, the total statistical uncertainty on the SBF is obtaine d as the square sum 
of the fitting uncertainty on Pq, and the default 25% uncertainty on P^es ( iTonry et al.lll990l ). 
The effect of the sky uncertainty is negligible on the SBF, while it is the main source of 
error in the color. In addition, these measurements suffer for a total systematic uncertainty 
due to: 1) the PSF normalization 0.03 mag), 2) the fit of LFs (~ 0.05 mag), 3) the filter 
zero points (~ 0.01 mag), and 4) if necessary, the transformation from the ACS photometric 
system to the standard system 0.02 mag). Summing in quadrature all the systematic 
sources of erro rs, we find that a to tal systematic uncertainty ~ 0.1 mag affects our SBF 
measurements (ICantiello et al.ll2005l ). On average, this corresponds to a systematic error of 
~ 6% on distances and on the Hubble constant. If we also add the maximum systematic 
uncertainty of ~ 0.25 mag due to the presence of shells, the total systematic error becomes 
~ 12%. It is worth emphasizing here, again, that the 0.25 mag systematic uncertainty due 
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to the shells is an upper limit, as it has been derived from the worst case, i.e. farthest distant 
galaxy, prominent shells. 

Table [2] reports the final SBF measurements and statistical uncertainties for all galaxies. 
In the Table we also report the SBF and colors obtained transforming the ACS filters F475W 
and F814W into standard B and I passbands. 



Discussion 



4.1. Calibration, SBF models and distances 



The application of SBF method as distance indicator requires the calibration of the 
absolute SBF magnitude, M, versus the broadband color. To date, no calibration is available 
for the ACS F814W bandpass versus {FAYbW — F814:W)q color, in the color range of the 
present galaxies sample. Moreover, even transforming the {FA75W — FSIAW)^ color to the 
sta ndard (B — /)n c olor, the empirical calibration of SBF in the standard I band determined 



by ICantiello et al.l (l2005l ) is defined in the color range 2.0 < {B — /)o < 2.25, while all our 
galaxies have {B — I)q < 2.0 mag. Thus, we decided to extend to bluer colors the calibration 
of absolute SBF magnitudes, by using models. 

The SBF models used here are based on the most updated version of the code SPoT 
(Stellar Population Tools, Raimondo et al. 2005l lR. These models have the advantage of 
fitting the SBF and color of eUipticals, as well a s of resolved and unresolved stellar clusters, 
for a wide range of ages and metallicities (see iRaimondo et al.l and references therein for 
details). For the specific purpose of this study, we compute d the theoretical SBFs in the 
ACS VEGAMAG photometric system using the BaSeL 3.1 JWestera et aPbooi IPatriceUi 
20061 ) stellar atmospheres library. The resulting SBF magnitudes are reported in Table [31 



We start with the empirical calibration by I Jensen et al.l (120031 ): 



Mi = (-1.58 ± 0.08) + (4.5 ± 0.25) x [{V - /)o - 1.15] 



(3) 



defined in the color range 0.95 < {V — I)o < 1.3. This relation has a high degree of 
reliability, since the zeropoint magi iitude has bee r i obt ained using the improved period- 
lur ninosity relat i ons fo r Cepheids by lUdalski et al.l (Il999l ). and the slope is the one derived 
by iTonry et al.l (119971 ) based on group membership of galaxies. On the theoretical side, we 



http: / /www. oa-teramo. inaf.it /SPoT 
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fit the SBF versus color models with a robust straight line fitting process that minimizes 
the mean absolute deviation. We selected stellar population models in the age range 1.5 < 
t{Gyr) < 14, and the metallicity range —0.4 < [Fe/H] < 0.3 (Fig. panel a). As a result 
we obtain: 



Mi = (-1.6 ± 0.1) + (4.5 ± 0.2) x [(V" - J)o - 1.15] 
in very good agreement with Eq. [31 



(4) 



The same models, but for the ACS photometri c system, are c ompared to the obser- 
vational data from the ACS V irgo Cluster Survey (IMei et al.l 120071 ). We have adopted a 
DM~31.1 for the Virgo cluster (IFerrarese et al.ll2000l ). and transformed models into the AB- 
MAG photometric system for sake of homogeneity with data (Fig. E]). A linear fit to models 
provides a slope a = 1.4 ±0.1 and intercept f3 = 29.0 ± 0.1, in good agreement with the 
empirical fits provided by Mei and collaborators, who give a = 1.3±0.1 and (3 = 29.09±0.04. 
We have also checked the case of a double linear fit, as suggested by the authors. We obtained 
a = 0.8 ± 0.1 and (3 = 29.0 ± 0.1 in the color interval 1.0 < {F475W - F850LP)o < 1.3; 
while a = 1.7 ± 0.1 and (3 = 29.0 ± 0.1 in the color range 1.3 < {F475W - F850LP)o < 1.6. 
In both cases the fit from models agrees within la with the empirical calibrations. 

Before going further, we evaluate the capability of the SBF method to derive accurate 
distances of shell ellipticals, as their peculiar morphology might distur b the SBF rn e asure- 
ment. To this purpose, we selected a sample of shell galaxies from the iTonry et al.l (120011 ) 
database, for which the distance modulus (DM) is available from methods not related to 
the SBF technique (Table H]). The absolute SBF magnitudes, Mnon-sbf dms, of these shell 
galaxies, calculated as the difference between the SBF apparent magnitude, rfiToi and the 
non-SBF DMs listed in Table HI are re ported in Fig. [ 71 Th e solid line in the figure represents 
the empirical calibration, Eq. [31 from I Jensen et al.l (120031 ). Taking into account the uncer- 
tainties, the observational data are consistent with the empirical calibration. The median 
difference between the absolute SBF magnitude predicted using Eq. [31 and Mnon-sbf dms 
is: (Mjo3 — MjyoN-SBF DMs) = —0.06 ±0.22 mag, so that no systematic offset is recognized. 
Moreover, the data-point nicely overlap with models. These results provide a further support 
in using the SBF method to derive distances of shell ellipticals. 

Re lying on these agree ments, and using the same set of stellar population models match- 
ing the [Jensen et al.l (120031 ) equation, we derived: 



Mf814W = (-0.94 ± 0.20) ± (2.2 ± 0.2) x [{FA75W - F814H^)o - 2.0] 
the relation is shown in Fig. [51 (panel b) . 



(5) 
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Following the same procedure, we derived the Mj versus (B-I)o calibrat ion (see Fig. [5l 
panel c), which extends to colors {B — I)o < 2.0 the previous calibration from lCantiello et al. 
(120051 ). The fit to models provides: 



Mi = (-1.6 ± 0.2) + (2.1 ± 0.2) [{B - J)o - 2.0] . 



(6) 



Using Eq. 013 and the data in Table [21 we obtained the DMs reported in Tabled (Col. 
2 and Col. 3). 

We note that the latter two relations may suffer of systematic uncertainties which typi- 
cally affect stellar population synthesis models as, for example, the adopted library of stellar 
atmosphere models, especially for cool and bright stars. In Fig. \5\ (panel d) we show the 
changes expected if a different stellar atmospheres library is used to obtain stellar population 
models in the F475 W and F81AW filters . The new atmospheres library is the combination 
of stellar models by IWestera et al.l (120021) for relatively hot s tars (Te > 4500 K), and by the 
PHOENIX models for cool stars (iBrott fc Hauschildtl l2005l ). Adopting these atmosphere 
models, Eq. [5] becomes: 

Mpsuw = (-0.86 ± 0.20) + (2.3 ± 0.2) x [{FA75W - F8UW)o - 2.0]. 

The distance moduli obtained using the latter equation agree within uncertainties with 
those reported in Tabled On average we find that the new DM varies less than 0.05 mag. 

After transforming the {B — I)o colors in Table [2] to {V — J)o color, we also derived the 
distance moduli of the four galaxies using the empirical calibration, E g. ([3]) and the mj data 
in the Table [2l To this purpose, similarly to ICantiello et al. J2005h . we derived the (V-I)o 
versus (B-I)o color transformation, using the same set of models adopted to obtain the Eqs. 
HB By fitting a straight line to the models, we derived: 

{V - I)o,transf = 0.47 ± 0.02 X (5 - /)o + (0.21 ± 0.02) 

The distance moduli obtained from this procedure are in Col. (4) of Table O All the 
DMs from the different calibrations are in good agreement within uncertainties. 

The uncertainties reported in Table [5] come from the propagation of statistical uncer- 
tainties already described in Sect. [3l and from the calibrations uncertainty. In addition, 
the following systematic uncertainties should also be taken into account: i) ~0.1 mag from 
fiat-fielding, filter zeropoint, etc., and, ii) the upper limit ~0.25 mag which is the maximum 
uncertainty possible due to possible residual shells. 



Finally, we compare our distances with the ones obtained using the Hubble law, reported 
in Table [U The recession velocity adopted for each galaxy are corrected for the Virgo 
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+ Great Attractor + Shapl ey Supercluster infall, based on the local velocity field model 



given in iMould et al.l (120001 ) . The good matching of our DMs with the kinematic DMs 
{{\DMsBF — -DM/^qI < 0.2) confirms that it is possible to measure reliable SBF magnitudes 
for galaxies up to 100 Mpc, even in the optical bandpasses, with an uncertainty of ±8 Mpc 
(statistical) ± 6 Mpc (systematic, ± 10 Mpc if the upper limit uncertainty due to possible 
unremoved shells is taken into account). 



4.2. Hn determination 



Since our galaxies are located beyond 4000 km/s~^, they constitute a good sample to 
determine the Hubble constant, Hq, given that the effect of local deviation from the smooth 
Hubble fiow is minimized at this redshift. Using the distances based on the calibration 
Eq. [5l we estimate i^o = 76 ± 6 (statistical) ±5 (systematic) [±8 systematic including the 
upper limit uncertainty from possible unsubtracted shells] km/s/Mpc. It should be noted 
here that the Hq value reported has been derived adopting the SBF measurements and the 
theoretical calibrations presented in this work, i.e. it does not suffer for the uncertainty of the 
Cepheids calibration. On the other hand, the adopted calibration suffers for the uncertainties 
and ass umptions that are tip ically embedded in stellar population synthesis models (see for 
example ICharlot et al.lll996l ). However, the reliability of the present mo dels is tested against 
known observational data (for example the empirical calibration from I Jensen et al.l ). Such 
comparisons suggest that the theoretical systematic uncertainties are not larger than the 
quoted uncertainties. 

Even if this determination is based on only four galaxies, it is interesting to note that 
the Hq value derived is in good agreement with the final result from the HST Key Project 
Team, Hq = 72 ± 4 (statistical) ±8 (systematic) km/s/Mpc, and with the value Hq = 
70 ± 5 (statisti cal) ±6 (systeraatic) km/s/Mpc obtained by the same authors using only 
SBF distances (iFreedman et al.ll200ll ). Finally, such value also agr ees with the recent values 
Hq = 73±4(statistical) ±5 (systematic) km/s/Mpc determined by lRiess et al.l (120051 ). using 
the multicolor light curve shape metho d on two SNIa, an d Hq = 72 ± 4 (statistical) ±4 
(systematic) km/s/Mpc determined by IWang et al.l (120061 ) . obtained by using a sample of 
109 SNIa and the ACu method. 
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5. Summary and conclusion 



We presented F81AW SBF measurements from ACS images of four distant sliell elliptical 
galaxies with radial velocities between 4000 and 8000 km/s. By using the SBF method, the 
distance moduli of these galaxies are derived for the first time. We provided new calibration 
relations of the absolute SBF magnitude versus the integrated color, specifically the Mpsuw 
versus {F4:75W — F814:W)o. Moreover, the Mr versus (B — I)q relation presented here 



extends to bluer colors the calibration of ICantiello et al.l (120051 ). The calibrations are based 



on new SBF models computed with the SPoT code for the ACS and standard filters. These 
models are aimed to simulate single-burst stellar populations of age ranging from t = 1.5 
Gyr up to t = 14 Gyr and metallicity from Z=0.008 to Z=0.04. The use of a theoretical 
calibration is important not only because it is free from the uncertainties affecting empirical 
secondary distance indicators, but also because no empirical calibration of SBF magnitudes 
for these photometric bands is available in literature. To verify the reliability of the Mpsi^w 
versus {F475W — F814:W)o calibration we used exactly the same set of models to derive 



Jensen et al. 


2003; 


Mei et al. 


2007) 



empirical calibrations shows an extremely good agreement. This result a dded to positive tests 



already performed on resolved and unresolved stellar populations (e.g. 


Brocato et al. 


1999. 


2000; 


Cantiello et al. 


2003; 


Raimondo et al. 


2005; 


Faeiolini et al. 


2007 


) make us confident 



that the theoretical calibration presented here can be safely adopted to derive distances. 

On the observational side, our measurements suffer for ~0.1 mag systematic uncertainty 
in the SBF measurements coming from the filter zeropoint, flat fielding, and PSF normal- 
ization. In addition, the presence of possible unsubtracted shells can affect the estimation 
of SBF amplitudes; we estimated an upper limit uncertainty of ~0.25 mag from the worst 
case in the present data (PGC6240: the farthest galaxy, prominent shell). 

Coupling the SBF measurements with the theoretical calibrations we find distances in 
agreement with the ones obtained using the Hubble law. The present measurements enlarge 
the sample of galaxies beyond 100 Mpc with optical SBF distances. 

Finally, using our SBF distances, we derive Hq = 76 ±6 (statistical) ±5 (systematic) [±8 
systematic when the upper limit uncertainty from possible unsubtracted shells is included] 
km/s/Mpc, in good agreement with the most recent estimations of the Hubble constant. 



It is a pleasure to acknowledge B. Patricelli for making available the color-temperature 
transformation for the ACS filters computed for her Laurea thesis. We also thank the ref- 
eree for her/his useful suggestions which improved the paper. Financial support for this 
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work was provided by PRIN-INAF 2006 "Prom Local to Cosmological Distances" (P.I.: G. 
Clementini) . 
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Fig. 1. — Panels a): F814W/F475W band surface brightness profile as function of r^/"^ (open 
circles/triangles). Panels b): measured (B-I)o profile as function of r^/^. Data are corrected 
for Galactic extinction and K-correction. 




Fig. 2. — The LFs of external sources. Open squares mark observational data; model fits 
to globular cluster and galaxy LFs are shown as dashed and dotted lines, respectively, and 
their sum as solid line. 
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Fig. 3. — Prom left to right: the original I-band frame, the residual frame and residuals 
times the final adopted mask, for each galaxy (upper left quote in each panel). 
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Fig. 4. — For each labeled galaxy, the top panels show the azimuthal average of the power 
spectrum (the average radii of the annuli used are reported in each panel) . The observational 

data (dots) are fitted by the sum of a scaled PSF power spectrum plus the constant white 
noise term (solid line). The lower panels show the Pq obtained as a function of the starting 
wavenumbcr of the fit, Po^kfu > k). The final Pq adopted is the weighted average of values 
in the fiattcr Pq vs k region. We excluded the lowest /c-numbcrs, that have been corrupted 
by the background subtraction, and the highest /c-numbers, that have been corrupted by the 
drizzling procedure. 
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^F475W-F814W)r 



Fig. 5. — Absolute SBF magnitudes versus integrated colors. In the figure we plot our 
models of different metallicities (as labeled) and ages (t=1.5, 2, 3, 4, 5, 7, 9, 11, 13, and 
14 Gyr). Symbols with increasing size mark models of older age. Full triangles mark the 
observational data. In each panel we report the linear fit to models (solid line) and the 
corre lation coeffici e nt (upper left label). Panel a: The calibration of Mi versus {V — /)o 
from lJensen et al.l (120031 ) is shown with a dot-dashed line. Panel b: As in panel (a) but the 
absolute Mpguw magnitude and {F475W — FSIAW)^ color are used. Panel c: As in panel 
(a) but the {B — I)q color is used. Panel d: As in panel (6) but a different stellar atmosphere 
library is used (upper right label). 
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(F475W-F850LP)o 



Fig. 6. — mpsbOLP versus {FA75W — F850LP)q. In the figure we plot the new SBF models in 
the ACS photometric system (ABMAG), shifted to a Virgo dista nce modul u s of 3 1.1(symbols 
are as in Fig. [5]). The observational data (full dots) are from iMei et al.l (120071 ). The best 



fit of the present SBF models is shown as a solid line, the best fit of the data by iMei et al. 
(120071 ) is reported as a dashed line. 
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Fig. 7. — Mj versus {V — Symbols are as in Fig. [5l In the figure we plot the cal- 
ibration from iJensen et al.l (120031 . solid line), la dispersion lin e s are reported with dashed 
line. A sample of shell elliptical galaxies from iTonry et al.l ( 120011 ) are reported as tri- 
angles. The absolute SBF magnitudes are derived coupling the apparent SBF measure- 
ments, with the SBF-independent DMs reported in Table HI It is noticeable the absence 
of any systematic deviation of shell ellipticals with respect to the empirical calibration 
((Mj03 — Mnon-sbf DMs) = —0.06 ± 0.22), and to the stellar population models. 
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Table 1: Properties of galaxies. 



Galaxy 


T 


R.A. 


Dccl 


^rec (.km/ s) 


DM 




Aj^S14W 


Af475W 


F814W 


F475W 




















exposure time 


exposure time 




















(sec) 


(sec) 


(1 ) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


PGC6510 


-3 


26.591 


-83.400 


4650 ± 22 


34. 0± 0.2 


0.588 


0.246 


0.492 


2712 


7874 


PGC 10922 


-2 


43.400 


-83.142 


4825 ± 42 


34.1 ± 0.2 


0.370 


0.156 


0.312 


3970 


9350 


PGC 42871 





191.022 


-34.202 


6400 ± 36 


34.7 ± 0.2 


0.290 


0.121 


0.242 


5908 


8369 


PGC 6240 


-3 


25.379 


-65.615 


7936''± ■ • ■ 


35.2 ± 0.5 


0.088 


0.036 


0.072 


7300 


21440 



NOTE- Col. (1): Galaxy name. Col. (2): Morphological T-type from RC3. Cols.(3-4): Right Ascension and 
Declination from RC3 (J2000). Col. (5): Recession velocity in the CMB reference frame corrected for Virgo 
+ Great Attractor + Shapley's infall (data from the NED database, www. http://nedwww.ipac. caltech.edu/). 
Col. (6): Kine matic Distance Modu lus obtained using Hq = 73±5 km/s/Mpc (from NED). Col (7) B-band 



extinction from 
using 



Sch 



Sirianni et al 



egel et all (|l99a i. Col. (8-9): F814W- and F475W-band extinctions, calculated from Ab 
(2005|) coefficients. Col. (10-11): Total exposure time for F814W and F475W images. 



Tor this galaxy only the Virgo infall corrected recession velocity is available. 
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Table 2: Color and SBF measurements. 



Galaxy 


Annulus 


Po 


Pres 


™/814W,0 


(F475W-F814W)o 


in/,0 


(B-I)o 




(arcsec) 


(ADU/s) 


(ADU/s) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


iA) 


(5) 


(6) 


(7) 


(8) 


PGC6510 


6-35 


0.0031 ± 0.0001 


0.0002 


31.65 ± 0.04 


1.47 ± 0.05 


31.70 ± 0.04 


1.86 ± 0.05 


PGC 10922 


13-40 


0.00225 ± 0.00004 


0.0002 


31.93 ± 0.04 


1.42 ± 0.04 


31.97 ± 0.04 


1.79 ± 0.04 


PGC 42871 


13-24 


0.0045 ± 0.0002 


0.0005 


32.18 ± 0.04 


1.35 ± 0.03 


32.22 ± 0.08 


1.70 ± 0.03 


PGC 6240 


16-30 


0.00125 ± 0.00006 


0.0001 


32.63 ± 0.06 


1.23 ± 0.03 


32.68 ± 0.06 


1.57 ± 0.03 



NOTE- The magnitudes are extinction and K-corrected. Col. (1): Galaxy name. Col. (2): Average 
annulus inner-outer radii for SBF measurements. Col. (3) Average Pq ■ Col. (4): Unmasked external sources 
fluctuation contribution P^e^,- Col. (5): SBF magnitude. Col. (6): (F475W - F814l^)o color. Col. (7): 
SBF magnitude in the standard photometric system. Col. (8): {B — I)q integrated color. 
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Table 3: SBF models. 



Age F435W F475W F55()M F555W F606W F625W F7T5W F814W F850LP (F475W - F814W)o 
(Gyr) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



Z = 0.008 



1.5 


1.988 


1.201 


-.021 


.277 


-.361 


-.766 


-2.075 


-2.387 


-3.090 


1.298 


2.0 


2.097 


1.287 


.086 


.376 


-.249 


-.642 


-1.959 


-2.277 


-2.987 


1.395 


3.0 


2.167 


1.363 


.200 


.479 


-.129 


-.510 


-1.871 


-2.203 


-2.928 


1.501 


4.0 


2.332 


1.532 


.382 


.658 


.055 


-.322 


-1.708 


-2.050 


-2.793 


1.504 


5.0 


2.242 


1.465 


.353 


.619 


.038 


-.328 


-1.645 


-1.973 


-2.696 


1.575 


7.0 


2.264 


1.475 


.368 


.632 


.060 


-.312 


-1.711 


-2.061 


-2.809 


1.621 


9.0 


2.429 


1.664 


.587 


.843 


.277 


-.079 


-1.450 


-1.805 


-2.569 


1.656 


11.0 


2.380 


1.636 


.587 


.835 


.284 


-.064 


-1.433 


-1.805 


-2.601 


1.694 


13.0 


2.431 


1.698 


.658 


.905 


.354 


.008 


-1.398 


-1.776 


-2.579 


1.72 


14.0 


2.393 


1.668 


.639 


.883 


.337 


-.007 


-1.412 


-1.795 


-2.605 


1.737 


Z = 0.01 


1.5 


2.094 


1.323 


.113 


.412 


-.230 


-.634 


-2.014 


-2.347 


-3.086 


1.309 


2.0 


2.270 


1.471 


.271 


.563 


-.066 


-.460 


-1.835 


-2.177 


-2.935 


1.42 


3.0 


2.402 


1.594 


.414 


.701 


.080 


-.305 


-1.732 


-2.094 


-2.879 


1.528 


4.0 


2.489 


1.690 


.527 


.810 


.196 


-.185 


-1.657 


-2.039 


-2.854 


1.573 


5.0 


2.304 


1.506 


.368 


.643 


.048 


-.322 


-1.712 


-2.074 


-2.862 


1.612 


7.0 


2.350 


1.572 


.468 


.734 


.154 


-.206 


-1.608 


-1.983 


-2.783 


1.648 


9.0 


2.470 


1.707 


.630 


.888 


.325 


-.029 


-1.394 


-1.771 


-2.584 


1.687 


11.0 


2.465 


1.701 


.626 


.883 


.317 


-.036 


-1.427 


-1.807 


-2.621 


1.725 


13.0 


2.477 


1.737 


.690 


.940 


.389 


.042 


-1.346 


-1.739 


-2.581 


1.763 


14.0 


2.572 


1.833 


.782 


1.034 


.477 


.129 


-1.289 


-1.694 


-2.553 


1.766 


Z = 0.02 


1.5 


2.579 


1.848 


.647 


.956 


.311 


-.091 


-1.582 


-1.994 


-2.903 


1.431 


2.0 


2.723 


1.933 


.708 


1.018 


.369 


-.030 


-1.478 


-1.889 


-2.816 


1.524 


3.0 


2.903 


2.096 


.874 


1.183 


.539 


.145 


-1.352 


-1.786 


-2.734 


1.615 


4.0 


2.926 


2.139 


.964 


1.260 


.642 


.262 


-1.185 


-1.628 


-2.606 


1.673 


5.0 


2.803 


2.006 


.832 


1.126 


.613 


.137 


-1.283 


-1.704 


-2.629 


1.703 


7.0 


2.912 


2.121 


.974 


1.259 


.662 


.294 


-1.074 


-1.489 


-2.415 


1.749 


9.0 


3.003 


2.229 


1.108 


1.385 


.800 


.438 


-.906 


-1.320 


-2.260 


1.793 


11.0 


2.981 


2.204 


1.080 


1.357 


.769 


.407 


-.910 


-1.324 


-2.271 


1.826 


13.0 


2.992 


2.230 


1.123 


1.396 


.814 


.455 


-.901 


-1.337 


-2.322 


1.862 


14.0 


2.979 


2.230 


1.147 


1.413 


.846 


.493 


-.820 


-1.248 


-2.227 


1.876 


Z = 0.04 


1.6 


2.903 


2.198 


1.036 


1.340 


.734 


.350 


-1.135 


-1.619 


-2.699 


1.633 


2.0 


3.057 


2.318 


1.146 


1.460 


.846 


.465 


-.999 


-1.482 


-2.564 


1.604 


3.0 


3.308 


2.646 


1.373 


1.674 


1.073 


.696 


-.744 


-1.248 


-2.390 


1.691 


4.0 


3.428 


2.669 


1.489 


1.789 


1.192 


.819 


-.614 


-1.121 


-2.277 


1.743 


5.0 


3.073 


2.303 


1.152 


1.444 


.863 


.499 


-.819 


-1.269 


-2.318 


1.785 


7.0 


3.245 


2.453 


1.295 


1.587 


1.005 


.643 


-.701 


-1.171 


-2.247 


1.846 


9.0 


3.273 


2.494 


1.360 


1.645 


1.073 


.717 


-.600 


-1.059 


-2.132 


1.886 


11.0 


3.315 


2.662 


1.442 


1.719 


1.160 


.809 


-.443 


-.890 


-1.968 


1.922 


13.0 


3.405 


2.660 


1.649 


1.824 


1.266 


.917 


-.330 


-.772 


-1.842 


1.956 


14.0 


3.281 


2.629 


1.438 


1.709 


1.157 


.811 


-.449 


-.902 


-1.990 


1.97 



NOTE- Col. (1): Age. Col (2-10): absolute SBF magnitudes in various ACS photometric filters. Col(ll): 
(F475W-F814W)o integrated color. 



- 28 - 



Table 4: Non-SBF DMs of shell galaxies and SBF magnitudes from iTonry et al.l (120011 ) 



Galaxy 


rriToi 


DM 


Distance Indicator 


Reference 






(mag) 


(mag) 








(1) 


(2) 


(3) 


(4) 


(5) 




NGC1316 


29.83 


31.2 ±0.1 


PNLF'^ 


Ferrarese et al. 


^2000) 


NGC 1344 


29.67 


31.4 ±0.2 


PNLF 


Teodorescu et al. 


(2005) 


NGC3923 


30.26 


31.5 ±0.2 


GCLF 


Silfkema et al. ( 


2007) 


NGC 4278 


29.34 


31.1 ±0.1 


GCLF 


Kundu & Whitmore r2001) 


NGC 5128 


26.05 


27.7 ±0.2 


Cepheids 


Ferrarese et al. 


^2007) 


NGC 4552 


29.39 


30.84 ±0.09 


GCLF 


Kundu k Whitmore r20011 



NOTE- Col. (1): Galaxy name. Col. (2): apparent SBF magnitude measured by iTonrv et all (|200lh . 
corrected for extinction. Col (3): DM obtained from the distance indicator listed in Col(3). Col (4): 
Reference for the DM. 
"Planetary Nebula LF. 
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Table 5: Distance Moduli. 



Galaxy 


DMacsi'' 


DMsTD,BI^ 


DMsTDJensen^ 


(1) 


(2) 


(3) 


(4) 


PGC 6510 


33.7 ± 0.25 


33.6±0.25 


33.6 ± 0.15 


PGC 10922 


34.2 ± 0.25 


34.1±0.25 


34.0 ± 0.15 


PGC 42871 


34.7 ± 0.25 


34.6±0.25 


34.5 ± 0.15 


PGC 6240 


35.2 ± 0.25 


35.2±0.25 


35.1 ± 0.2 



NOTE- Col. (1): Galaxy name. Col. (2): DM obtained with eq. H Col. (3): DM obtained with eq. El 
Col(4): DM obtained with eq. [3] 
"Statistical uncertainties are reported. The systematic error is ^ 0.1 mag (see text). 



